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SUI>fliAEY 

The results of an investigation of \vater-vapor condensation 
shocks in the air passing through supersonic tunnels are presented. 
Criterions for condensation-free flow are established by correlating 
experimental observations with the Volmer theory of nuclei forma- 
tion. Experimental observations were made at Mach mmbers up to 
2.01. The criterions are presented in a form independent of 
tunnel-inlet stagnation pressure and are extended theoretically to 
a Mach number of 4.00. Preliminary evidence of the effect of tun- 
nel size on the criterions is presented. 


INTRODUCTION 

When a mixture of water vapor and air is expanded beyond sat- 
uration conditions in passing through a supersonic wind tunnel, a 
change of phase of the vapor resulting in a nonisentropic discon- 
tinuity of the flow is often observed. The discontinuity 1s 
referred to as a "condensation shock." Condensation shock, as used 
herein, however, refers only to that region of the air stream in 
which the increase in static pressure accompanying phase transition 
of the water vapor is measurable . 

Tlie change in static-pressure distribution, if it occvirs in 
the nozzle of a supersonic tunnel, is indicative of a deviation 
from the ordered expansion of the flow. The design of the ensuing 
part of the nozzle is thereby invalidated and nonparallel flow 
exists at the nozzle outlet. Nonparallel-flow affects, accompanying 
condensation shocks, have been proven to complicate the interpreta- 
tion of aerodynamic data (reference l). 

A theory for the formation of condensation nuclei that appears 
to be applicable to the treatment of condensation shocks is pre- 
sented by Tolmer in reference 2. This theory is modified herein 
and used in con.junctlon with experimental investigations made at 
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the MCA Lewis laboratory in order to establish criterions for 
condensation-free flow. The criterions are presented in a form 
independent of tvinnel- inlet stagnation pressure. The physical 
phenomena related to the condensation are briefly discussed. 


ANALYSIS 


The partial pressure of a vapor p^ expanded through a super- 
sonic tiuinel may be assigned to change in accordance with the adia- 
ba.tic gas law. The satmat ion- vapor pressure p^, however, may 


be asstmied to change in accordance with the Clausius-Clapeyron 
relation. During expansion of the air through a timnel, the 
partial pressure of the vapor' therefore decreases as a power func- 
tion of the temperature, whereas the saturation- vapor pressure 
decreases as an exponential function of the temperatvire. Thus, 
with adiabatic expansion, the partial pressure of the vapor may 
exceed the saturation- vapor pressure; that is, Py/^ ^ which 


results in a supersaturation of vapor for a given temperature. 
This condition may be relieved by the transformation of the 
excess vapor into the liquid phase as droplets. 


A requisite for the formation of droplets is nuclei upon 
which vapor can condense. A common source of nuclei may be the 
hygroscopic substances in the atmosphere. In a supersonic tunnel, 
however, the amount of vapor condensing on such available nuclei 
during their rapid passage through the tunnel was calculated to be 
negligible by Oswatitsch (reference 3). Condensation of water 
vapor in the absence of foreign nuclei has been observed with the 
Wilson cloud chamber (reference 4). In this reference it has been 
postulated that aggregates of water molecules form the required 
condensation nuclei. Because heavy condensation shocks have often 
been observed in supersonic tunnels, nuclei of this type may be 
assumed to exist at, or Just prior to, the condensation shock. 


In a system where aggregates of water molecules are contin- 
ually forming and breaking up, an aggregate continually grows if 
the vapor pressure at the surface of the aggregate is less than the 
partial pressure of the vapor in the surrounding atmosphere. At 
equilibrium the radius of the droplet is defined in terms of the 
partial pressure of the vapor and the saturation- vapor pressure by 
the Thomson relation (reference 2, p. 19) : 

^ ^ m 1 

Pj Et logg (p^/d^) 
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where 


r radius of aggregate, centimeters 

S surface tension of liquid, dynes /centimeters 

p density of liquid, grams/cubic centimeter 

V 

m molecular weight of vapor, grams/raole 

E universal gas constant, joules/(®C) (mole) 


t absolute temperature, °K 

p /i^ ratio of partial pressin:*e of vapor to saturation- vapor pres- 
sure in equilibrium with infinite liquid surface 


A theory for determining the rate of nuclei fomatlon is 
advanced by Volmer (reference 2). It is assumed that aggregates 
of the size defined by equation (l) are the condensation nuclei of 
the supersonic stream. Volmer 's equation is modified in appendix A 
to the following more convenient expression: 


logg J = logg 9.589 X 10^^ + 


logg 



k/A H 

p 


1/6 

~2jE 

I 


-»• 


2 logg p^/p^ 


17.558 



(loge 


( 2 ) 


where 

J rate of formation of nuclei, number per second per cubic 

centimeter 





ergs 


(All symbols used herein are defined in appendix A.) 


The estimated variation of Kg with temperature is discussed 
in appendix B and presented in figure 1. The saturation- vapor 
pressures p for supercooled water were determined from a relation 
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by Osborne and Meyers (see reference 5). The density of water, 
used to the 2/3 power in the second term of the right side of equa- 
tion (2), was obtained by extrapolation of the measured values 
given in reference 6, Because of the slight change in density with 
temperature, the effect of a.ny error in the density (used to the 
2/3 power) on the value of logg J was considered negligible. The 
variation of logg J of equation (2) with vapor-pressin*e ratio 
(expressed as logg P^/p^) for selected isotherms is presented in 

figure 2. The sensitivity of the nuclei formation rate to the 
vapor-pressixre ratio is characterized by the steep slopes of the 
isotherms (fig. 2); for example, for a gas temperature of -54.4° F, 
a change in nuclei formation rate from 1 to 20,000 per second 
per cubic centimeter corresponded to a change in vapor-pressure 
ratio from 10.3 to 13.2 

In order to associate the ‘relation of figure 2 with the prop- 
erties of a gas flovring through a supersonic tunnel, ciurves AB 
and BC are presented. If the temperature of the air entering a 
supersonic tunnel is increased, the vapor pressures p^ in the 

tunnel remain constant but the saturation- vapor pressures p 

00 

change. The resulting relation between local stream temperature t 
and the quantity logg (p^/p^) at a given local Ifech number is 

illustrated by curve AB (directed from A to B). The rate of 
nuclei formation at the station, and hence the tendency for con- 
densation, is decreased. If, having reached the conditions cor- 
responding to B for the given station, the humidity and the tem- 
perature of the inlet air are maintained constant and flow condi- 
tions at higher Mach numbers farther along the tunnel are considered, 
the curve BC is indicative of the accompanying trend in the ten- 
dency for condensation, the value of J increasing with an Increase 
in Mach number. Thus, if the value of J required for condensation 
is attained in the tunnel, increasing the temperature lowers the 
J value at each station and expansion to a higher I4ach number is 
required befoi-e the critical J is again reached. The experi- 
mental problem is to determine J for condensation-free flow. The 
experimental data in figure 2 are presented later. 


APPARATUS 


The tunnel system (fig. 
The specific humidity of the 
specific humidity of 0.00025 


3 ) is a nonreturn steady-flow type, 
inlet air can be varied from a minimum 
pound per pound of dry air to a specific 
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humidity corresponding to saturation. The inlet-air dew points 
(hereinafter called dew points) were determined to within ±0.5° F 
hy simultaneous use of a continuously Indicating electronic dew- 
point meter (reference 7), a manually operated dew-point meter, and 
an expans ion- type dew-point meter. The inlet-air stagnation tem- 
perature (hereinafter called inlet temperature) can he regulated 
from 40° to 250° F hy means of thermostatic controls and was meas- 
ured with a 5/l6-inch aspirating thermocouple (reference 8) and a 
self -ha lancing potentiometer to within ±0.5° F. 

Tunnels with 3.4- hy 3.4-inch and 4- hy 10-lnch test sections 
were used in the investigations. Three interchangeable nozzles 
were available for the 3.4- hy 3.4-incb tunnel. One nozzle vras 
designed for a 14ach number of 2.00 with a straight diverging wall 
(included angle of 3.50°). The throat was faired into the plane 
diverging- wall section by use of an arbitrary radius. This nozzle 
installation is shown in figure 4(a). The other two nozzles were 
designed hy the soi-irce-f low method (reference 9) for Mach numbers 
of 2.01 and 1.45 (figs. 4(h) and 4(c), respectively). The nozzle 
used in the 4- hy 10-inch tunnel was also designed by the source- 
flow method for a i4ach number of 2.00 (fig. 4(d)). An approximate 
correction for viscosity effects was made in the design of the 
source-flow nozzles in order to attain the design Mach number in 
practice. 

Each tiinnel was eq.uipped with wall orifices of approximately 
0.013-inch diameter throughout the expanding surfaces and side 
plates. Typical side-plate instrumentation is presented in fig- 
ure 4(e). The orifices, as well as the pitot rakes located at the 
inlet and the outlet of each tunnel, were coupled with multiple- 
tube mercury manometers that could be simultaneously photographed. 
The change in height of any mercinry column could be read to the 
nearest 0.005 inch by means of special viewers. The orifices of 
the 4- by 10-inch tunnel could also be coupled with a multiple- 
tube oil manometer, on which a change in column height correspond- 
ing to 0.0006 inch of mercury could be detected. 


BESULTS AMD DISCUSSION 
Condensation-Free Criterions 

The increase in static pressure of a supersonic stream 
passing through a condensation shock is theoretically predicted in 
appendix C and is experimentally established in reference 1. This 
increase in static pressure can be used to determine the location 
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of a condensation shock in a supersonic tunnel. An example of the 
experimental data used in determining the shock location is shown 
in figvire 5. The variation in static-pressure ratio p/Pq with 
inlet temperature Tq at one arbitrary station in a supersonic noz- 
zle is presented for a dew point of 56° F. (Dew point is defined 
as the temperature at which the liquid, or solid, and vapor phases 
coexist for a given partial pressure of the vapor. Temperatirres 
for the remainder of the report will he given in English units.) 

The shock was moved along the nozzle in the direction of increasing 
Mach number by raising the inlet temperature. At an inlet temper- 
ature of 160° F (point A in fig. 5), a minimum static-pressure 
ratio of 0.201 was reached and all measurable effects of the con- 
densation shock were downstream of the station. An inlet tempera- 
ture of 160° F is the minimum required for condensation-free flow 
at the given station and dew point and is designated the critical 
temperature T^. 

If the critical temperature corresponding to each station 
throughout the supersonic nozzle and the test section of a tunnel 
is established, the position of the condensation shock as a function 
of the inlet temperature is determined. Critical temperatures were 
obtained at stations throughout the nozzles of the 3.4- by 3. 4- inch 
tunnel by varying the inlet temperature and maintaining the dew 
point constant at 18° to 23.5°, 49°, and 56° ±1° F (table I, inves- 
tigations 1, 2, and 3, respectively). These critical temperatures 
are presented in figure 6 as a function of the experimental Mach 
numbers measured at the tunnel stations. The conditions required 
for condensation-free flow at a given dew point must fall below 
the curve coiresponding to the dew point. The shock movement 
throughout the test section of the tunnel of 1.45 Mach number 
(3.4- by 3.4-in. test section) is presented in figure 7 in a manner 
similar to that for the nozzles of figure 6 for dew points of 10.8°, 
20.0°, 26.8°, and 44.5° F (table I, investigations 4, 5, 6, and 7, 
respectively) . In recording the movement of the shock through the 
test section, however, the location of each station (referenced 
from the test-section inlet) is used (fig. 7) in place of Mach num- 
ber because the Mach number throughout the test section is nearly 
constant. The shock was moved down the nozzle (in the direction 
of increasing Mach number) and through a considerable part of the 
constant-area test section by increasing the inlet temperature. 

An inlet temperature was ultimately attained at which the shock 
appeared to move suddenly out of the test section; that is, the 
critical temperatures for the stations dotnistream of the last 
observed shock position were reached simultaneously. This inlet 
temperature Tg (fig. 7) corresponds to the data that appear to 
approach vertical lines. 
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The stream temperature and the vapor-pressure ratio and thus 
the nuclei formation rate corresponding to condensation-free flow 
in the test section of the 5.4- by 5.4-inch tunnel were established 
for each inlet stagnation temperature Tg of figure 7 and are 
presented in figure 2. The stream temperature was obtained from 
Tg and the energy relation 

where M is the ^•^ach number at the test-section inlet. The vapor 
pressure of water p^ at the stream tempei^ture was calculated by 

means of reference 5. The partial pressure of the vapor p^ at the 

test-section inlet was determined by using a combination of Dalton's 
law of partial pressures and the perfect gas law. The specific humid- 
ity a in pounds of water vapor per pound of dry air is therefore 


a = 0.6223 — (4) 

P - Py 

Measiired values of inlet pressure and dew point were used to obtain 
the inlet specific humidity. Any change in the specific humidity 
prior to the condensation shock was considered negligible. The 
stream pressure p iras measured with wall orifices. 

The experimental procedure used in investigations 4 to 7 
(table l) in the 3.4- by 3.4-inch tunnel was also used for inves- 
tigations 8 and 9, which were made at a Mach number of 2.00 in the 
4- by 10- inch tunnel. 

The critical temperatures throvighout the constant-area test 
section are presented in figure 8. The curves in figure 8 do not 
appear to approach vertical asymptotes as do those of figure 1 , 
but are inclined with a positive and nearly constant slope. This 
positive slope is more in keeping with the nuclei theory of con- 
densation than is the infinite slope of figure If because as the 
stream temperatiare increases and the rate of formation of conden- 
sation nuclei decreases it is reasonable to believe that a greater 
time (and thus a greater distance) would be required for the com- 
pletion of the condensation process. An explanation of the appar- 
ently infinite slope of figure 7 may be that the slope of test- 
section distance for a given change in inlet temperature at the 
high dew points of investigations 4 to 7 is too great to be 


8 


MCA TR 2518 


detected in a tunnel having the limited length of the 3.4- hy 

3. 4- inch tunnel. Also, the boundary- layer growth, known to 
have a greater effect on the flow properties in small tunnels 
than in large tunnels, may have restrained the condensation in 
the test section through an increase in the saturation pressure 
of the vapor accompanying the deceleration of the flow. 

If the first 14 Inches of test section of the 4- by 10- inch 
tunnel is considered to be the test section, and the nuclei forma- 
tion rates corresponding to the critical temperatures at 14 Inches 
from the test-section inlet (fig. 8) are calculated in the same 
way as ^ras done for the inlet temperature Tg in the 3.4- by 

3.4- inch tunnel, the average nuclei formation rate compares well 
with the average rate required for condensation-free flow in the 

3.4- by 3. 4-inch tunnel (fig. 2). 

A condensation-free-criterions chart based on the VoLmer 
theory for nuclei formation may be constructed from figure 2. The 
chart will apply to tunnels comparable in size to the 3.4- by 

3.4- lnch tvinnel and also to tunnels comparable in size to the 4- 
by 10- inch tiumel and having test sections approximately 14 inches 
in length. The average of the six experimentally determined values 
of logg J (fig. 2) was selected as the criterion for condensation- 
free flow and designated logg Jg. The variation of stream temper- 
ature with vapor-pressure ratio required for condensation-free flow 
(fig. 9) corresponding to this criterion was obtained from a cross 
plot of figure 2 at logg J = logg Jg. 

• 

The relation between the properties of the tunnel air and the 
criterion may be illustrated with the curve AB of figure 2, as 
presented in figure 9. If the gas temperature at a station in the 
tunnel is increased (moving from A to B) , the corresponding 
decrease in the degree of saturation is such that the curve passes 
to the left of the criterion curve into the region of condensation- 
free flow. 

The relation given in figure 9 may be expressed in terms of 
the more convenient tunnel operating parameters of temperature, Mach 
nvmiber, and dew point. Calculations were made for various values of 
streeim temperature t and inlet temperat\n*e T. The Mach number 
corresponding to any selected t/T was calculated using equa- 
tion (3). For a given value of t, logg (p^/p„j) was obtained from 

figure 9 and the saturation- vapor pressure p was evaluated in 

00 

accordance with the Osbome-Meyer relation (reference 5) . Hence 
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the partial pressure of the vapor Py was determined. The assumption^ 
Justifiable in moat cases, is made that py and Py are negligible 

compared with p and P, respectively, and that a is constant. 

It follows from equation (4) that Py/Py = p/P and thus the partial 
pressiare of the vapor at stagnation may be calculated from the 
isentropic relation 


A condensation-free-criterlons chart relating the partial 
pressure of the vapor at the tunnel inlet (expressed as dew point), 
the inlet temperatvire, and the test-section Mach number for the 
tunnels used in the investigation reported is presented in fig- 
ure 10, The inlet temperatures are indicated with solid lines and 
the stream temperatures used in the calculation are shown as 
dashed lines. Although the properties of supercooled water were 
used in the calculations, frost-point rather than dew-point tem- 
peratures were used on the chart below 32° P, because the usual 
dew-point meters read frost point in this range. 

The results of the experimental investigations used in deter- 
mining logg Jg are also presented in figure 10. The numbers by 
the symbols indicate the measured inlet stagnation temperature. 

The deviation of the experimental values of logg J from logg Jg 
in figure 2 correspond in figure 10 to dew-point fluctuations of 
approximately ±1.5° F or to inlet-temperature changes of approxi- 
mately ±2° F. 

Wo assumption was made regarding tunnel- inlet pressure in the 
calculation of flgin*e 10, except that the dew point be measured 
at the tunnel- inlet pressure. The partial pressure of the vapor 
(expressed as dew point), the specific humidity, and the mixture 
pressure (water-vapor pressure plus air pressure) are presented in 
convenient form in figure 11 in accordance with equation (4). 

Frost points rather than dew points were again used below 32° F. 
Figure 11 combined with the condensation-free-criterions chart 
(fig. 10) can be used to predict the tunnel- inlet-air properties 
required for condensation-free flow in the tunnel test section for 
tunnels comparable in size with the 3.4- by 3.4-inch and 4- by 
10- inch tunnels and operating over a range of inlet pressures. 

For example, consider a tunnel operating with an inlet pressure 
of 100 pounds per square inch absolute and having a test-section 



. -JL. 

, r-1 


(5) 
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Mach munher of 2.00. Assume that, because of the available drying 
apparatus, the minimum specific humidity of the air at 100 pounds 
per square inch absolute is 0.0002 pound per poixnd of dry air. 

From figure 11, a specific humidity of 0.0002 corresponds to a dew 
point of 10.5° F at an inlet pressure of 100 pounds per square 
inch absolute. From figure l5, the inlet stagnation temper*ature 
(solid lines) corresponding to a Mach number of 2.00 and a dew 
point of 10.5° F is approximately 228° F. This temperature cor- 
responds to the minimum value required for condensation-free flow 
in the tunnel test section. In practice, it would be advisable to 
exceed the minimum temperature as a protection against the possi- 
bility of local regions of higher Mach number "triggering" the 
conderxsation process. 

The condensation-free criterions of figure 10 were based on 
investigations made at a Mach nimiber of 1.45 in a 3.4- by 3.4-inch 
tunnel and at a Mach number of 2.00 in a 4- by 10- inch tunnel with 
a 14-lnch test section. An additional investigation (10, table l) 
at I«Iach number 2.01 and dew point 17.3° F was made in the 3.4- by 
3.4-inch tunnel to ascertain the effect of tunnel Mach nismber on 
the conditions necessary for condensation-free flow. An inlet- 
stagnation temperature of 246° F was required for condensation-free 
flow throughout the test section as compared with the inlet stagna- 
tion temperatvire predicted in figure 10 of 243° F. Although one 
investigation may not be considered conclusive, it does indicate 
that theoretical extension of the data for a Mach number of 1.45 
on a basis of nuclei-f oimiation theory is valid up to a Mach number 
of 2.01. Preliminary unpublished data of tunnels operated at sev- 
eral Mach nimibers up to and including a Mach number of 2.00 and 

having sizes li to 2 times that of the 4- by 10-inch tunnel, sub- 

La 

stantiate (±10° F inlet stagnation temperature for a given dew 
point and Mach number) the chart of figure 10. 


Scale Effect 

The extension of figure 10 to tunnels having test sections of 
considerable length, however, requires further consideration 
because of the scale effect apparent in figure 8. It was found 
that an increase in test-section length necessitated an increase 
in inlet temperature to avoid condensation. The Increase in inlet 
temperature per foot increase in test-section length for the 4- by 
10- inch tunnel (Mach number of 2.00) was 5.3° and 6° F for inlet 
dew points of -5° and -10° F, respectively. This obseorved scale 
effect agrees with preliminary unpublished NACA data for a tunnel 
5 times the size of the 4- by 10- inch tunnel. 


1089 


1089 


NACA TN 2518 


11 


A scale effect was also noticed In the difference between the 
movement of the shock throu^ the nozzles of 2.01 Mach nimiher of 
the 3.4- by 3.4-inch tunnel (dew point 17.3° F) and 2.00 Mach num- 
ber of the 4- by 10- inch tvumel (dew point 19° F), which may be 
attributed to the different rates of expansion of the flow in the 
two nozzles and the time req.uired for the condensation process. 

For a given temperature (fig. 12), the shock occurred at a higher 
^fech number in the nozzle of the 3.4- by 3.4-inch tunnel than in 
the nozzle of the 4- by 10-lnch tunnel, the 3.4- by 3.4-inch tun- 
nel nozzle expanding the flow to a 14ach number of 2.00 in less than 
one-third the distance required by the 4- by 10-inch tunnel nozzle. 
Hence the data obtained throiighout the nozzles are not generally 
acceptable for use in constructing criterion charts. 

The criterion logg was determined using vapor pressures 
over supercooled water rather than over ice, even though the 
stream temperatures in the region of the condensation shock were 
well below freezing. Although no conclusive evidence regarding 
the physical state of the condensation droplets has yet been 
obtained, it has been observed that ice formations on tunnel walls 
and models occur only under conditions of extreme moisture content 
and recovery temperatures close to freezing. An ice formation on 
the glass sides of the 3.4- by 3.4-inch tunnel (nozzle Mach num- 
ber, 1.45) is shown in figiire 13 for a dew point of approximately 
70° F and an inlet temperature such that, with 90-percent recovery 
of the kinetic energy, the wall temperature was a few degrees 
above freezing. 

Physical Appearance of Condensation Shock 

Schlieren photographs of condensation shocks are presented in 
flgiure 14. Figures 14(a) and 14(b) are photographs of a shock 
located in a channel having no curvature Just downstream of the 
throat, and figures 14(c) and 14(d) are photographs of a shock in 
a nozzle having considerable curvature from the region of the 
throat to the test section. In figure 14(a) the shock, located 
Just do\mstream of the throat, appeared well defined and extended 
over a relatively small distance; whereas in figure 14(b) the shock, 
having been moved into a region of higher velocity, appeared poorly 
defined and extended over a considerably greater distance in the 
direction c£ the flow. The shock in both positions appeared to be 
normal to the stream. In the nozzle having considerable curvature, 
the shock located Just downstream of the throat was also normal to 
the stream and well defined, as presented in figure 14(c). As the 
shock was moved downstream, it again required a greater distance 
to attain completion (fig. 14(d)) and at the same time \raa changed 
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from a shock normal to the stream to one having an X- shape. In 
figure 14 (d)^ well-defined disturbances appeared to be passed back 
and forth across the tunnel downstream of the shock. 


Practical Considerations for Tunnel Design 


The obvious way to insure flow that is free from condensation 
effects is to eliminate the condensation shock by drying or heating 
the air or both in accordance with criterions such as those pre- 
sented in figure 10. The requirements, however, become rapidly 
prohibitive as the tunnel Mach member increases. Specifications 
must therefore be established for the maximum allowable vapor con- 
tent of the inlet air at and below which condensation effects are 
negligible. No attempt has been made to determine precisely such 
criterions. Information was gathered during the investigations 
discussed herein, however, that may be used for a qualitative 
approach to the problem. The magnitude of the change in static 
pressure along the walls of the test- section inlets was observed 
for several specific humidities at Mach numbers of 1.45 and 2.01 
in the 3.4- by 3. 4 -inch tunnel and at a Mach number of 2.00 in the 
4- by 10-inch tunnel. The change in static pressure across the 
shock for the two Mach numbers is presented in figure 15 as a 
function of the specific humidity. The difference between the 
average pressure ratio behind the shock (pressure ratio at B in 
fig. 5) and the pressure ratio corresponding to the minimum 
pressure ahead of the shock (pressure ratio at A in fig. 5) for 
stations at the test-section inlet was averaged to arrive at the 

P2 - Pi 

change in static-pressiire ratio — 5 of figure 15. The 

Po - Pi ^ 

theoretical — was calculated for Mach numbers of 1.45 and 

Pq 

2.00 as a function of specific humidity in accordance with appen- 
dix B. It appears from experiment and theory that for a given 


P9 “ Pi 

liach number, increases with an increase in specific 

Po 


humidity and that the magnitude of 
humidity decreases with an increase 


^2 “ ^1 

— for a given specific 

^0 

in test-section Mach number. 
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SUMMAEY OF RESULTS 

From an experimental investigation of condensation shock in 
supersonic tunnels, the following results were obtained: 

1. The Volmer theory of nuclei formation and an experimental 
determination of the rate of formation of nuclei required for a 
condensation shock were combined to establish criterions for 
condensation-free flow in a 3.4- by 3. 4- inch tunnel and a 4- by 
10- inch tunnel with a 14- inch test section. 

2 . An increase in inlet temperature was required to avoid 
condensation when the test section of the 4- by 10-inch tunnel 
was lengthened, inlet dew point and Mach number remaining con- 
stant. This scale effect was found to agree with unpublished 
NACA data for a tunnel five times the size of the 4- by 10-inch 
tunnel. 

3. The static-pressxire rise across a condensation shock 
located in the test- section inlet of a tunnel increased with an 
increase in specific humidity for a given test-section Mach 
number, and, for a given specific humidity, decreased with an 
increase in Mach number, as predicted by theory. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, June 16, 1949 
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APEENDIX A 
SYI-IBOLS 

The following symbols are used in this report: 
specific heat at constant pressure 
specific heat at constant volume 
rate of formation of nuclei, number /( sec ) (cc) 



Mach number 

molecular weight of vapor, gram /mole 
stagnation pressvire 

partial pressure of vapor at stagnation 

static pressure 

partial pressure of vapor 

saturation- vapor pressure in equilibrium with infinite liquid 
surface, mm Hg 

ratio of partial pressure of vapor to saturation- vapor pres- 
sure (relative humidity) 

heat added to flow at condensation shock, ft-lb/lb dry air 

universal gas constant, joules /(°C) (mole) 

radius of aggregate, cm 

surface tension of liquid, dynes/cm 

inlet-air stagnation temperatvire 

critical temperature (inlet-air stagnation temperature corres- 
ponding to condensation-free flow at a given station) 
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Tg Inlet-alr stagnation temperatvire corresponding to 
condensation-free flow in tuimel test section 

t stream temperatvire 

V Telocity, ft/sec 

y ratio of specific heats, Cp/Cy 

p density, gr€uns/cc 

p density of liquid, grams /cc 

L 

Q specific humidity, Ih water vapor/lh dry air 

T duration of state of greatest super saturation, sec 

Subscripts : 

0 before tunnel inlet 

1 before condensation shock 

2 after condensation shock 
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APPEMDIX B 


MODIFICATION OF VOLMER THEORY 

A theory for determining the rate of nuclei formation is pre- 
sented by Volmer in reference 2. The probability of the fonaation 
of nuclei is evaluated therein (p. 136) as 


logg Jt = logg 


T + logg (8.5 X 10^^) + logg 



+ 



(Bl) 


where 
9.5 X 10^^ 

17.49 


combination of physical constants, 
combination of physical constants, 


deg 

_ (dyne) (cm)_ 

(mole^) (deg^) 
(dyne^) (cc) 


The combinations of the physical constants were corrected 
herein to agree with the currently accepted values of the universal 
gas constant and Boltzmann's constant. Accordingly, 9.5 x 10^^ ■vras 
changed to 9.589 X 10^^ and 17.49 was changed to 17.558. 

Durii]g the exiansion of the air throu^ the tunnel and the 
accompanying change in stream temperature and vapor- pres sure 
ratio, the only significant change in the value of the terms on 
the right side of equation (Bl) occinrs in the last term, now cor- 
rected to 
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The inlierent diff j.culty in applying the Volmer theory is apparent 
upon examination of this term when it is realized that the correct 
values of siurface tension and density of water must he determined 
for stream temperatures well below freezing. It is proposed that 
surface tension, density, and molecular weight be combined as 


for purposes of extrapolating their knora values, inasmuch as the 
rate of change of Kg with temperature is nearly constant over the 
kno^m range. A plot of Kj; calculated from the physical properties 
of water in accordance with reference 6 is presented in figure 1. 
Because of the difficulty of measuring the properties of super- 
cooled water (the admitted possibility for error is shown by the 
shaded area), the experimental upward trend of Kg below 273° K 
was neglected and the data above the freezing point were extra- 
polated linearly to 170° K. If Kj; is inti’oduced into equa- 
tion (B1), and the duration of critical conditions logg t is 
dropped from both sides of equation (Bl), the eqiiation may be 
Witten 



logg J = loge 9.589 X 10^^ + logg 




(B2) 
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APPENDIX C 


CliANGE IN STATE ACROSS CONDENSATION SHOCK 


The lisual flow relations may be used to obtain a one-dimensional 
approximation of the change in state across a condensation shock. 

The following assumptions greatly simplify the calculations (see 
following figure) : 


The shock takes place at constant area, the specific heats before 
and after condensation do not change, the gas constant of the air is 
not appreciably affected by the moisture of the air, and all the 
moisture condenses out at the shock until saturation is reached. 

The energy relation may be written 



Condensation shock 



Throat 



^ 2 ^ 

+ q, - Cp tg + -Y“ 


(Cl) 


Relation (Cl) may be rearranged into the form 


tg 7Eti " 7-1 





(C2) 


1 

7 - 1 2 



2 


With the momentum relation. 


Pi + PlV = P2 + P2^2^ 


(C3) 
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and the perfect gas law, 


p = pEt 


(C4) 


and the assumption of conservation of mass, relation (C2) may he 
rearranged to give 


1 MgS V 

7^ T~ _ y - ^ 2 

V~% / V~% / 


The value of may he determined from the criterions chart 
(fig. 10) if the specific humidity and the inlet temperature are 
known. The correct value of q, is obtained hy a series of approx- 
imations (temperatvire and q are interdependent). Hence M 2 is 
determined. The effect of the condensation shock upon the stream 
pressure may then he determined from equation (C3), reanranged to 
give 

P, 1 + 7WL.2 

^ \ (C6) 

Pq 1 + 7M2^ 


Equation (C6) and the isentropic relation 




(C7) 


may he used to calculate the change in static pressure as presented 
in figure 15. 
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TABLE I - INVESTIGATIONS OF CONDENSATION 
IIJ SUTERSONIC FLOW 


Investi- 

gation 

Test-section 
Mach number 

Dew point 
(OF) 

Tvinnel 

(in.) 

1 

No test section 

18-23.5 

3.4 by 3.4 

2 

2.01 

49 

3.4 by 3.4 

3 

No test section 

56ii 

3.4 by 3.4 

4 

1.45 

10.8 

3.4 by 3.4 

5 

1.45 

20.0 

3.4 by 3.4 

6 

1.45 

26.8 

3.4 by 3.4 

7 

1.45 

44.5 

3.4 by 3.4 

8 

2.00 

-5 

4 by 10 

9 

2.00 

-10 

4 by 10 

10 

2.01 

17.3 

3.4 by 3.4 

11 

2.00 

19 

4 by 10 



170 190 210 230 250 270 290 310 330 

Absolute ten5)erature , t, 


Figure 1. - Variation of Kg with ten^erature. 
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Vapor -pres sure ratio, logg(p^/p^) 

Figure 2. - Variation of rate of formation of nuclei with Tap or -pres sure ratio for Tarlous 

Isotherms* 


Rate of formation of nuclei, logg J 
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4-21-48 



id) 4- by 10- inch tunnel installation show- 
ing nozzle; Mach number, 2,00. 



C-21290 

4-21-48 


(a) 


(b) 


3.4- by 3.4-inch 
tunnel nozzle In- 
stallation with 
included angle of 
3.50°; nozzle in 
place; Mach num- 
ber, 2.00. 



C-21254 

4-21-48 



3.4- by 3.4-inch 
tunnel nozzle; 
Mach number, 
2 . 01 . 



C-21284 

4-21-48 



(e) 4- by 10- inch tunnel Installation show- 
ing instrumentation. 


c/zT^l 

4-21-48 


Figure 4. - Installations and nozzles of 3.4- by 3,4-lnch and 4- by 10- inch tunnels 




Static-pressxire ratio, p/P 
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Figure 5. - Variation of static-pressure ratio with inlet temperature for given station. 

Dew point, 56° P. 
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Mach number 


1089 



F*lgupe 6« - Relation between Mach, mmber and critical temperatures throughout nozzles of 

3«4- by 3#4-inch tunnel at three dew points# 
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Dlsteuice from test-section inlet, in 
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• - Critical ten?)eratures throu^out test section of 3.4- by 
3.4-incli tunnel. Macb number, 1.45. 


Figure 7 
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Critical inlet tempera txire , T^, ®F 


Figure 8. - Critical temperatures throu^out test section of 4- by 
10-inch tunnel, Mach number » 2.00. 
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0 1 2 3 4 5 


Vapor-pressxare ratio, 

• - Variation of stream temperat-ure with vapor -pres sure 
ratio. Logg J, 4.1. 


Figure 9 
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Test-section-inlet 



Figure 10 


- Criterions for 
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condensation-free flow 



MJjcture pressure, in, Hg 


34 


MCA IW 2518 



NACA TN 2518 


35 


Mixture 



water vapor/ lb dry air 


specific humidity for constant dew points. 



Ivlach ntunber 


CKl 

cn 


Tunnel Macsli 
(in*) number 


O 3.4 by 3.4 2.01 

O 4 by 10 2.00 



Critical inlet temperature, 


Figure 12. - Critical temperatures for various Mach numbers throughout nozzles of 
3.4- by 3.4-lnch and 4- by 10-inch tumels at Mach numbers of 2.01 and 2.00, 
respectively. Dew point, 17. 3^ to 19. 0^ p. 
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Figure 13. - Ice formation In 3.4- 'by 3.4-inch tunnel. Mach number, 1.45; dew point, 70° F. 





C-21438 
5- 14-48 


C-21436 
5- 14-48 


(a) Condensation shock at throat in slovly diverging 

nozzle. 


(b) Condensation shock in figure 14(a) showing change 
in thickness with movement downstream. 




C-21439 
5- 14- 48 


(c) Condensation shock at throat in rapidly diverging 

nozzle. 


(d) Condensation shock in figure 14(c) approaching 
conventional X- shape with movement downstream. 


Figure 14. - Schlleren photographs of condensation shocks in supersonic flow in nozzles of 3.4- by 3.4-lnch tunnel (flow from 

left to right). 
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NACA-Langley - 12 - 6-51 - 1000 



Figure 15. - Change in static-pressure ratio with specific humidity for condensation at 

Mach numbers 1, 45 and 2, 00 ^ 
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